The smallest possible heat engines 
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We construct the smallest possible self contained heat engines; one composed of only two qubits, 
the other of only a single qutrit. The engines are self-contained as they do not require external sources 
of work and/ or control. They are able to produce work which is used to continuously lift a weight. 
Despite the dimension of the engine being small, it is still able to operate at the Carnot efficiency. 



Very recently considerable progress has been made in 
understanding foundational aspects of thermodynamics 
by addressing a new class of questions: whether there 
exist additional fundamental limitations on thermal ma- 
chines, arising specifically due to their size [1]; whether 
there is an absolute minimum size to all thermal ma- 
chines; alternatively, if machines are small, is their per- 
formance constrained, for example can they still achieve 
the Carnot efficiency? 

For the case of refrigerators it was shown [1] that there 
is no minimum size when self-contained machines con- 
sisting of only two qubits, or a single qutrit, were dis- 
covered. Furthermore it was shown that these refrigera- 
tors have two important properties; they are able to cool 
towards absolute zero [1J and they can operate at the 
Carnot efficiency |2|, showing that here size appears not 
to limit the machines at all. 

One advantage in focusing initially on refrigerators 
rather than heat engines is because they allowed us to 
avoid the explicit notion of work in our study, whilst 
making progress in understanding small thermal ma- 
chines. Work is a concept that appears difficult to cap- 
ture in the quantum regime, where ideas such as "or- 
der" usually associated to it are hard to make sense of. 
Therefore it is not clear that the existence of small self- 
contained refrigerators implies anything about small 
self-contained heat engines. 

The study of the thermodynamics of quantum sys- 
tems has become a thriving and dynamic field HHU- 
The definition of work has indeed been discussed exten- 
sively [7-9] and applied to quantum heat engines IfTOl — 
16] and Carnot cycles fT7ti22l . The focus however has 
been on quantum thermal machines which explicitly or 
implicitly have macroscopic objects in the background 
which supply either work or some form of control; for 
example systems which are externally driven, or make 
use of sequences of unitary evolutions. The definitions 
of work given so far have been applicable and very well 
suited to this situation. Our focus is however differ- 
ent; we are interested here in self contained heat engines, 
where there is no external work or control - the only 
external interaction being with thermal reservoirs. The 
previous definitions of work therefore do not apply di- 
rectly to our situation, and an alternative must be found. 



In this work we will use exactly the definition of work 
put forward by Carnot, namely: 

"motive power (work) is the useful effect that a 
motor is capable of producing. This effect can al- 
ways be likened to the elevation of a weight to a 
certain height. " fl3j 

to construct a quantum heat engine. Our engine will 
produce work in the same way as above, by lifting a 
weight; here the weight consists of a system with an in- 
finite number of energy levels, and the work created by 
our heat engine causes the position of the weight to in- 
crease with time. We will present two models for the 
heat engine, one consisting of two qubits (two level sys- 
tems), and a second consisting of a single qutrit (three 
level system). We will also show that they can operate 
at the Carnot efficiency. 

TWO-QUBIT MODEL 

The heat engine consists of two qubits. Qubit 1 is in 
contact with a "cold" heat bath at temperature T c and 
qubit 2 is in contact with a "hot" bath at temperature 
T/j. Like every heat engine, this engine acts by extract- 
ing heat from the hot bath and converting it into work, 
while dumping some heat into the cold bath. 

We imagine that the engine delivers work by pulling 
up a weight. The weight is isolated from both baths. 
To simplify the situation we consider that the weight is 
pulled up very slowly, so that we can neglect the change 
of its kinetic energy and consider only the potential en- 
ergy. Furthermore, we suppose that the weight can be 
situated only at some discrete equidistant heights, so 
that the energy difference between them is the same. 
Hence, the weight has discrete energy eigenstates \n) w 
with corresponding energy eigenvalues E n = n£, with 
£ > 0. Alternatively, we can imagine that the engine 
delivers work by pumping energy into a harmonic os- 
cillator; both situations are formally almost equivalent 
(the harmonic oscillator energies are limited from below 
byO). 

In the absence of interaction we take qubit i, i = 1, 2 to 
have energy eigenstates |0), and and corresponding 
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FIG. 1: Schematic diagram of two-qubit heat engine. Qubits 
1 and 2, with energy level separations £j and E2 are in contact 
with thermal reservoirs at temperatures T c and 7), respectively. 
To this a weight is connected, with separation £ = E2 — E\. 
The particles interact with each other and the temperatures are 
chosen such that the transition where the weight is lifted (solid 
green arrows) is biased over the transition where the weight 
falls (dashed orange arrows). 



energy eigenvalues and E; > 0. The free Hamiltonian 
of the system is thus 

00 

H = Ei|1)i(1| + E 2 |1>2<1|+ E n£\n) w (n\ (1) 
The energies are taken such that 



E 2 - Ei = £. 



(2) 



Given this constraint, the energy levels |01, n) and 
1 10, n + 1) are degenerate. The engine acts by making 
transitions between these degenerate states. 

The qubits and weight interact via the Hamiltonian 

H in t=gE (|01,n)(10,n + l| + |10,n + l)(01,n|J. (3) 

n=— 00 

The basic idea behind the design of this engine is to 
bias the transition 1 01, n) — > \10,n + 1} in which the 
weight is lifted in favor of the reverse transition in which 
the weight is lowered. This is obtained by coupling 
the two qubits to heat baths at different temperatures, 
Th > T c , chosen such that the probability for the qubits 
to be initially in the state |01) is larger than the probabil- 
ity to be in the state 1 10) . 

In the absence of interaction each of the two qubits 
is at thermal equilibrium. Let the probabilities for the 
ground and excited state of qubit i be r, and rj. At ther- 
mal equilibrium they are related by 

T\ — r\e kJ c , r 2 — r 2 e kT h, (4) 

where k is Boltzmann's constant. This leads to the prob- 
abilities qio of the state 1 10) and (joi of |01) to be 



lio = V2 = nr 2 e 



(5) 



Therefore in the absence of interaction, the probability 
of the state |01) is larger than that of 1 10), i.e. (joi > <7l0' 
when we have 



(6) 



To demonstrate the functioning of our engine we also 
need to model the interaction of the two qubits with 
their heat baths. We take here a simple model, the same 
as the one used in [1J. Namely, we consider that in a 
short time interval St with probability pjSt the state of 
qubit i is reset to its thermal state, the density matrix t,, 



i i = r i |0) i (0|+?7|l)i(l|. 
The overall equation of motion is 



(7) 



do 2 
1 - -i [H + H int/ + E Pi (TfriP ~ P) (8) 

t=l 



dt 



where p denotes the state (density matrix) of the full sys- 
tem comprised of qubits 1 and 2 and the weight. This 
Master equation is consistent in the limit where g and 
Pi are small, otherwise there are corrections to the dissi- 
pative dynamics which must be take into account ||23] . 
In this regime we expect our device will act as a heat 
engine and continuously lift the weight as long as the 
relations Q and |6) are fulfilled. 

The quantity of interest is the average energy of the 
weight, (E w ). We expect that, after a transient period, 
this average energy will increase uniformly, that is we 
expect that 



dt 



const > 



(9) 



We now show that this is indeed the case. From <|Sj we 
obtain 



dt 



dt 



Tr(H w p) = -ig£A(t) 



(10) 



where 



A(0 = E(< 01 ' M li°l 10 ' M + l)-(10,n + l|p|01,n>) (11) 

n 

If we define further the two quantities T\(t) and T 2 (t) 
which are the instantaneous ground state probabilities 
for qubits 1 and 2 respectively, 



ri(t) = E(< 00 ' M li°l 00 ' n } + ( 01 ' M H 01 ' n }) (12) 

n 

r 2 (f) = J^({00,n\p\00,n) + {W,n\p\10,n)) (13) 

n 

then together, along with A(f), these three quantities 



3 



obey the coupled set of equations 



dt^ 



2ig(T 1 (t)-T 2 (t))-(p l + p 2 )A(t) 



dt 



dt 



Ti(t) = +igA(t) + p 1 (r 1 -T 1 (t)) 
T 2 (t) = -igA{t) + p 2 (r 2 -T 2 {t)) 



(14) 



This set of equations can easily be solved. The solution 
for A(f) is given by 



A(f) = X> 



i=l 



kit + ngpmin-n) 

(p 1 + p 2 )(2g 2 + p 1 p 2 ) 



(15) 



where the Sj are constants determined by the initial con- 
ditions. The eigenvalues A, are the solutions of the 
characteristic equation (A + Pi)(A + p 2 ) (A + p\ + p 2 ) + 
2g 2 {2\ + p\ + p 2 ) = 0, and all have negative real part. 



The first term in | |T5| | is therefore the transient behaviour 
which decays in time, and thus in the long time limit 
only the second term survives and A becomes a con- 
stant. Thus, from equation (TlOj, 



lim — < 

t-¥oo dt 



2£g 2 p 1 p 2 (r 1 - r 2 ) 
{p 1 +p 2 )(2g 2 + p 1 p 2 ) 



(16) 



which is positive whenever r\ > r 2 , which is equivalent 
to cjiQ > qoi ~ the condition (|6) imposed initially to en- 
sure a probability bias in the correct direction. Therefore 
the average energy of the weight increases linearly with 
time, thus demonstrating that the weight is able to be 
continuously lifted. 

We note however that the probability distribution of 
the position of the weight, (and therefore the probabil- 
ity distribution of its energy) will be expected to also 
spread in time. Essentially we expect the evolution, af- 
ter the transient period, to be that of a biased random 
walk. This can be confirmed to be the actual behaviour 
of the system. A similar (but more involved) calculation 
to the above shows that AE 2 , goes like t in the limit that t 
tends to infinity. Thus the standard deviation in the en- 
ergy, AE W , varies asymptotically as as expected for a 
random walk. 



THE EFFICIENCY 

The most important property of the engine is the ef- 
ficiency with which it can run, namely the amount of 
work that can be delivered for a given amount of heat 
extracted from the hot bath. 

The rate at which heat flows between the qubits and 
their environments is given by the change in energy of 
each qubit due to the interaction with the baths. Given 
the Master equation |8|, the heat currents are therefore 



dt 



Qi^PiTrlHiixi-piit)) 



(17) 



Recall that Ti(t) and T 2 (t) are the ground state popu- 
lations of p\{t) and p 2 (t) respectively. The asymptotic 
solution for these two variables is given by 

limr,-(Q = PlP2(Pi+P2fo+y(Pin + P2r 2 ) (18) 
(Pi + P2){2g 2 + Pip 2 ) 

Substitution of this result into \17) , and making use of 
|TJ and it is found that 

±Q. = (-l)i 2E iZ 2 PtP2(n-r 2 ) ( 9) 

dt Ul [ lj ( P i + p 2 )(2g2- [U) 



indicating, as expected, that heat is flowing out of qubit 
1 into the 'cold' reservoir and is flowing into qubit 2 
from the 'hot' reservoir. Identifying now 4r (E w ) as the 
work extracted from the engine, since this is the energy 
which is used to raise the weight, it is seen that the effi- 
ciency rp- of the quantum heat engine is given by 
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dt ' 



E 2 



dtQi 



E 2 



E 2 



(20) 



However, for the weight to be lifted, condition Q must 
be satisfied. Thus, in the limit where the speed at which 
the weight is lifted becomes infinitely slow, that is when 
||6} becomes an equality, the quantum fridge reaches its 
maximum possible efficiency fjmaxr given by 



7n 



(21) 



which is exactly the Carnot efficiency of a heat engine 
running between two reservoirs at temperatures T c and 
T/,. Thus no new additional constraints arise due to the 
small size of the engine presented, with the only con- 
straint upon its efficiency being the Carnot limit. 



ONE-QUTRIT MODEL 

There is a yet smaller self contained heat engine than 
the one described above. It consists of a single qutrit 
- a quantum system with a three dimensional Hilbert 
space. Imagine that the three transitions (from |0) to |1), 
from |0) to 1 2), and from |1) to |2)) are affected by dif- 
ferent heat baths. This might occur, for example, due 
to the physical geometry of the system comprising the 
qutrit. Given such a set-up, it is possible to construct a 
smaller heat engine. The overall set-up can be seen in 
the Schematic diagram in Fig. [2] 

The corresponding Master equation for this system 
can also be written down and solved for the long-time 
behaviour of the system. The detailed calculation can 
be found in the Appendix. The solution, in the situation 
where all thermalisation rates are equal, denoted here 
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FIG. 2: Schematic diagram of qutrit heat engine. The qutrit 
is assumed to have each of its transitions at different tempera- 
tures. The temperature of the transitions is shown in the inset 
on the left. A weight is again connected, with energy level sep- 
aration £2, matching that of the separation of the two upper 
levels of the qutrit. The particles interact with each other and 
the temperatures are chosen such that the transition where the 
weight is lifted (solid green arrows) is biased over the transi- 
tion where the weight falls (dashed orange arrows). 
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by p, is found to be 



t^co at 



ig 2 E 2 p(r r r h -r c r h ) 



2g 2 (2 + r c + r r ) + 3p(3 - r c r r - r c r h - r r r h ) 



(22) 



where r c , r r and 77, are the ground state thermal popu- 
lations of qubits at the corresponding temperatures, de- 
fined as in Q, and g is again the interaction strength. 
The condition that the average energy of the weight 
asymptotically increases in time is now given by 



r r r h > r c r h 



(23) 



which says that, in the thermal states, the product of 
probabilities for state |2) must exceed that of the product 
of probabilities of state |1). 

Acknowledgements We acknowledge support from 
EU integrated project Q-ESSENCE. 



[1] N. Linden, S. Popescu and P. Skrzypczyk, Phys. Rev. Lett. 

105,130401 (2010). 
[2] P. Skrzypczyk, N. Brunner, N. Linden and S. Popescu, 

arXiv:quant-ph/1009.0865 
[3] S. Carnot, Reflexions sur la puissance motrice du feu (1824); 
Reflection on the Motive Power of Fire, Dover (1960). 
G. Gemma, M. Michel and G. Mahler, Quantum Thermo- 
dynamics, Springer (2004). 
[5] E. P. Gyftopoulos and G. P. Beretta, Thermodynamics: Foun- 
dations and Applications, Dover (2005). 
A. E. Allahverdyan, R. Ballian and Th. M. Nieuwen- 
huizen, J. Mod. Opt. 51, 2703 (2004). 



[4] 



[6] 



APPENDIX: ONE-QUTRIT MODEL DETAILS 

In the absence of interaction, the free Hamiltonian of 
the qutrit and weight is given by 

H = Ei|l)i(l| + (Ei +E 2 )|2)i(2| +E 2 J2n\n) a (n\ 

(24) 

The qutrit and weight interact with each other via the 
interaction Hamiltonian 

(25) 



H mt = g£(|l,n + l)(2,n| + \2,n){l,n + l 



The qutrit, due to its physical geometry, has each of its 
transitions in contact with a different thermal reservoir. 
The transition between the states |0) and |1) is at tem- 
perature T c , the transition between the states |0) and |2) 
is at the temperature T r and the transition between the 
states |1) and |2) is at the temperature T n . The Master 
equation governing the dynamics of the qutrit interact- 
ing with the environments and weight is given by 



dp 
dt 



-i[H + H int , P\ 
7 c (t c «0| P |0) + 



(26) 



where 



Jr«0|p|0) 

p h (r h ((l\p\l) 



r c \i)(i\, 



;i| j0 |l)) + |2)(2|(2|p|2)- |0 
;2| j0 |2)) + |l)(l|(l| j0 |l)-p 
(2|p|2)) + |0)(0|(0|p|0)-p 



r C |0)(0| 



T r = ry|0)(0| 
fh\2)(2\ 



fr\2)(2\, 
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are the thermal states for each of the three transitions 
respectively, where the r, are defined according to Q. 

To show that the qutrit engine is also able to produce 
work, the quantity of interest again is the asymptotic be- 
haviour of the mean energy of the weight. From ( |26) it 
follows that 



dt 



dt 



Tx{H w p) = -igE 2 n{t) 



(27) 



where 



n(f) = £ ({2,n\p\l,n + 1) - (l,n + l\p\2,n)) (28) 



then together with Q(f) they satisfy the following set of 
equations 



dt 
d 

dt 
d 

di 



n = +2ig(B 2 -B l )-(p c + p r + p h )Ci 



(30) 



Bi = -igCl + p c (r c (l-B 2 )-B 1 )-p h (rj l B 1 -r h B 2 ) 



B 2 = +igCi + p r {r r {l-B 1 )-B 2 ) 



{r h Bi - r h B 2 ) 



Similarly to the two-qubit case, by defining further the 
two variables B\ (t) and B 2 (t), the instantaneous excited 
state probabilities of the qutrit, 

B x {t)=Y d {l,n\p\%n), B 2 (t) = J2(2,n\p\2,n) (29) 



The behaviour is found to be qualitatively the same as 
the two-qubit heat engine; after an initial transient pe- 
riod, Q(f) tends to a constant. Thus the long-time be- 
haviour of the mean energy of the weight is given by 



lim -j-AEw) 
t-yoo dt 



2g 2 E 2 (p c p r (r r - r c ) - p h (r h - r h )(p c r c + p r r r )) 



2 g 2 (Pc{l + r c ) + Pr{l + r r )) + {p c + p r + Ph)(PcPr{l ~ f~ c 7>) + p c Ph{^ ~ T^h) + PrPh(^ ~ r r r h )) 



(31) 



and thus the weight is lifted whenever the following 
condition is satisfied 



h ( r r - r c ) > p h {r h - r h ) ( p c r c + p r r r ) . (32) 



